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ABSTRACT. The existence of bright and dark multi-bump solitary waves for Ginzburg- 
Landau type perturbations of the cubic-quintic Schrodinger equation is considered. The 
waves in question are not perturbations of known analytic solitary waves, but instead arise 
as a bifurcation from a heteroclinic cycle in a three-dimensional ODE phase space. Using 
geometric singular perturbation techniques, regions in parameter space for which 1-bump 
bright and dark solitary waves will bifurcate are identified. The existence of TV-bump dark 
solitary waves (N > 1) is shown via an application of the Exchange Lemma with Exponen- 
tially Small Error. TV-bump bright solitary waves are shown to exist as a consequence of 
the work of Kapitula and Maier-Paape [?] . 

1 Introduction 

The nonlinear cubic-quintic Schrodinger equation is given by 

iA t = A xx + \A\ 2 A + a\A\ i A, (1.1) 

where A is a complex- valued function of the variables (x,t) £ R x R + . When a = 0, 
the equation becomes the focusing cubic nonlinear Schrodinger equation, and is used to 
describe the propagation of the envelope of a light pulse in an optical fiber which has a 
Kerr-type nonlinear refractive index. For short pulses and high input peak pulse power 
the refractive index cannot be described by a Kerr-type nonlinearity, as the index is then 
influenced by higher-order nonlinearities. In materials with high nonlinear coefficients, such 
as semiconductors, semiconductor-doped glasses, and organic polymers, the saturation of 
the nonlinear refractive-index change is no longer neglible at moderately high intensities 
and should be taken into account ([?]). Equation ( |L~l| ) is the correct model to describe the 
propagation of the envelope of a light pulse in dispersive materials with either a saturable 
or higher-order refraction index ([?], [?]). 

Equation Ql.lD cannot really be thought of as a small perturbation of the cubic nonlinear 
Schrodinger equation, as it has been shown that a physically realistic value for the parameter 
a is \a\ ~ 0.1 ([?]). It turns out that the most physically interesting behavior occurs when 
the nonlinearity is saturating, i.e., a < 0, so for the rest of this paper it will be assumed 
that a ~ —0.1 ([?], [?], [?], [?], [?]). An optical fiber which satisfies this condition can be 
constructed, for example, by doping with two appropriate materials ([?]). 

One of the more physically interesting phenomena associated with the double-doped 
optical fiber is the existence of bright solitary wave solutions (|^4(x)| — * as \x\ — * oo) 
in which the peak amplitude becomes a two-valued function of the pulse duration. These 
solutions have two different peak powers, and were proven to be stable ([?], [?], [?], [?]). 
The solution with the lower peak power corresponds to a perturbation of the one that exists 
for the cubic nonlinear Schrodinger equation, while the one with the larger peak power is 
due to the saturating nonlinearity. 

Equation ( |1.1| ) describes an idealized fiber; therefore, it is natural to consider the per- 
turbative PDE 

iA t = (1 + ieai)A xx + ieaA + (1 + iedi)\A\ 2 A + {a + ied 2 )\A\*A, (1.2) 

where 0<6<1 and the other parameters are real and of 0(1). The parameter a\ describes 
spectral filtering, a describes the linear gain or loss due to the fiber, and d\ and d 2 describe 
the nonlinear gain or loss due to the fiber. So that ( |l.2|) is a well-defined PDE for e > 0, it 
will henceforth be assumed that a± > 0. 



1 



Solitary wave solutions to are found by setting 

A(x,t)=A(x)e ifMt , (1.3) 
and then finding homoclinic solutions for the ODE 

(1 + ieai)A" + (fj, + ieo)A + (1 + ied x )\A\ 2 A + (a + ied 2 )\A\ /L A = 0, (1.4) 
where ' = d/dx. Multiplying the above by 1 — iea\ and setting 

x = (1 + e 2 a\)x 

yields the equivalent ODE 



A" + ((p, + eVcr) + ie(cr - am)) A + ((1 + e 2 aidi) + ie(d 1 - ai))|v4| 2 ^4 
+((a + e 2 a 1 d 2 ) + ie(d 2 - aia))\A\ 4 A = 0, 



(1-5) 



where now ' = d/dx. 

To simplify matters, all the 0(e 2 ) terms in the above equation will be dropped. This 
step is taken only so that the notational complexity is made as simple as possible, and can 
clearly be done without any loss of generality, as for e small these terms are negligible. 
Upon dropping these terms, the ODE that will be studied can finally be written as 

A" + (/x + ie(a - am)) A + (1 + ie{d x - a x ))\A\ 2 A + (a + ie(d 2 - aia))|A| 4 A = 0. (1.6) 



Equation fll.q ) has been extensively studied by many authors ([?], [?], [?], [?], [?], [?], 
[?]). These papers have been concerned with finding various types of solutions, including 
fronts (kinks), bright solitary waves, and dark solitary waves. The methods employed have 
been both geometric ([?], [?], [?], [?]) and analytic ([?], [?], [?]). 

Bright solitary waves exist when there are solutions to (|0]) which are homoclinic to 
\A\ = 0. When e = with 

i^ii? (L7) 

one can find a bright solitary wave with a peak power larger than that associated with the 
cubic nonlinear Schrodinger equation. When e = and \i is in the range 

</x<0 (1.8) 
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there exists dark solitary wave solutions. The dark solitary waves are solutions to ( |1.6[ ) 
which have the property that |^4(x)| — > A§ ^ as |x| — > oo. 

It can be hypothesized that for the perturbative PDE fll.2j ) there may be a competition 
between the bright and dark solitary waves for \x sufficiently close to the critical value 

"•-is- (L9) 

It may be further hypothesized that for /j, sufficiently near /j,* it may be possible to construct 
novel solutions by gluing together the bright and dark solitary waves in some way. The goal 
of this paper is to explore this possibility. 

Doelman [?] has recently considered the problem of finding A^-circuit solutions. These 
TV-circuits are constructed by piecing together N copies of a dark solitary wave, and exist 
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as a solution to (1.6) for e sufficiently small and the parameters in a certain domain in 



parameter space. Recent work by De Bouard [?] has shown that the dark solitary waves 



are an unstable solution to (LI), due to the presence of a real positive eigenvalue for the 
operator obtained by linearizing about the wave. By a result of Alexander and Jones [?], 
since the original solitary wave has an unstable eigenvalue it is expected that the iV-circuit 
solution will have at least N unstable eigenvalues. 

It was previously stated that the bright solitary wave is a stable solution to ( |1 . 1[ ) . 
However, recent numerical work by Soto-Crespo et al [?] shows that this wave becomes an 
unstable solution to (|1.2| ) for e nonzero. The numerics suggest that this instability arises 
from the presence of a real eigenvalue for the linearized problem bifurcating out of the origin 
and into the right-half of the complex plane. 



Due to these results, when attempting to construct stable solitary waves to (1.2) one 



would hope to avoid using either the bright or dark waves which exist for (LI). This may 
by possible by considering the situation when // = //*. When e = and \x is equal to this 
critical parameter, there exists a front solution (kink) to ( |l.6[) which has the asymptotic 
behavior 

JO, x — > — oo 



\A(x) 



x — ► oo. 



Furthermore, since ( |L6j ) is invariant under x — ► —x, there also exists a solution which 
satisfies 

J Aq, x — > — oo 
I 0, x — > oo. 



\A(x) 



Thus, in the ODE phase space there exists a heteroclinic loop, from which solitary waves 
may bifurcate as e is made nonzero. By following the proof in De Bouard [?] it can be con- 
jectured that for the kink solution the linearized operator possesses no unstable eigenvalues. 
Therefore, it may be possible that any solitary waves bifurcating out of the heteroclinic loop 
may also not have any unstable eigenvalues. 

Based upon the above discussion, it will be of interest to study the ODE ( |L6| ) when // is 
(at least) 0(e) close to //*. In particular, it will be of interest to study the dynamics of the 
ODE near the heteroclinic cycle for e nonzero. Towards this end, it will first be shown that 
the cycle persists for e nonzero for fi = //(e), with //(e) — > //* as e — > 0. After then fixing 
e the dynamics will be studied as // bifurcates from //(e). It is important to note that // is 
a natural bifurcation parameter in this problem, as it is a free parameter (recall equation 

CD). 

This paper will be devoted to proving the following theorems. Set 

A °~ 4a' 

and define 

d\ = d\ — a\ 

d,2 = d-2 — OMi 

a = -(a + A^di + A$d 2 )/e. 

Assume that a = 0(1) for all e > 0. 

Before the main theorems can be stated, a preliminary lemma is needed regarding the 
persistence of the kink solitary wave for e / 0. 

Lemma 1.1 For < e -C 1 there exists a //(e), with |//(e) — //*| = 0(e 2 ), such that a kink 
solitary wave exists. 
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Remark 1.2 Since ( \1.6\ ) is invariant under the transformation x — > —x, the above lemma 
guarantees the existence of a heteroclinic cycle in the ODE phase space for e nonzero. 



Theorem 1.3 Let < e C 1, and assume that 

3 

(di + -A 2 d 2 )d < 0. 



There exists a /^(e) < A*( e )j with //(e) — /U/j(e) = 0(e C//<E ), suc/i i/iai a bright solitary wave 
exists. Furthermore, for each N > 2 there exists a bi-infinite sequence {^} such that when 



\i = there is an N -pulse solution to (1.6). The N -pulse is even in x. In addition, 



\^- fIh (e)\=0(e-^) 



as k — > oo. 



Remark 1.4 An N -pulse solution is constructed by piecing together N copies of a bright 
solitary wave. 

Remark 1.5 The work of Soto-Crespo et al /?/ suggests that the 1-pulse solution is indeed 



stable as a solution to (t.i) 



Theorem 1.6 Suppose that < f < 1, and that < \i — fi(e) = 0(e n ) for some n > 3. 
Further assume that the parameters satisfy 

(di + \Ald 2 )d < 
(di + 2Ald 2 )5 < 
{d l + Ald 2 ){di+pd 2 ) > 0, 

where 

P = Al + \Al(-ln^ + Otf) 

for 0<e<)j<l. Then there exists an N(e) > 1, with N(e) — > oo as e — > 0, suc/i i/ia^ 
N -circuit solutions exist for 1 < N < N(e). 

Remark 1.7 T/ie suitable domain in parameter space is shown in Figure^. In this figure 
the tildes have been dropped, [5 = a 2 ja\, and xq = Aq. 

The rest of the paper is organized as follows. In Section 2 the relevant three-dimensional 
ODE is set up, and the known solution structure is discussed. In Section 3 the relevant 
manifolds in the phase space are described. The bright, dark, and kink solitary waves 
are formed by intersecting the appropriate manifolds. This section is the most technical, 
and can be skipped on a first reading. Sections 4 and 5 discuss the intersection of these 
manifolds. In Section 6 the final proofs of the main theorems are given. 

Acknowledgements: I would like to thank Alejandro Aceves for introducing me to the 
cubic-quintic Schrodinger equation and some of the references pertaining to it. 
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2 Setup of equations 

stable bright solitary waves are possible for /j, near //*. They do not consider the case of 
dark solitary waves in their paper, however. 

In this section the first order system of ODEs that will be studied will be derived. The 
goal will be to write the system as compactly as possible so that the notational burden of 
the succeeding sections will be minimized. Define 



where 



and furthermore set 



A(r) = n* + r 2 + or 4 
oj(t) = d\r 2 + d?,r , 

d\ = d\ — a\ 
d2 = d2 — a\a, 

a = —{a — aifi) 
ji = n* -fl 



(2.1) 

(2.2) 
(2.3) 



(fi* is defined in (|Q|)). After setting 

A(x) =r(x)e i f^ ds 
and substituting this expression into (|1.6|) one gets the third order system 



r' = ru 



u' = -u 2 - A(r) + 4? + fi (2.4) 
(j)' = -2u4> + e(uj(r) - a). 

In the above equation the variable u is defined by u = r' jr. 

The above equations will now be simplified notationally. First, redefine the variables 

(x, r, u, 4>) by 

(t,x,y,z) = (x,r,u,(j)). (2.5) 
Next, remove the tildes from all the variables. The ODE ( |2.4| ) then becomes 

x' = xy 

y> = —y 2 + Z 2 — X(x) + H (2.6) 

z' = —2yz + e(to(x) — a), 

where now ' = d/dt and u>(x) = d\x 2 + c^x 4 . For the subsequent analysis it will be useful 
to set 

a = uj{xq) + ecr, 

where xq is the largest positive root of \{x) = (xq = — 3/4ct). Upon substituting the above 
expression into (2.6) and dropping the tilde, one obtains the equation 

x' = xy 

y> = _ y 2 + z 2 _ A ( x ) + ^ (2.7) 

z' = — 2yz + e(uj(x) — uj(xq) — ecr). 

Equation (^^) is the one which will be studied in the subsequent sections. However, it is 
important to notice that the equation has the symmetry given in the following proposition. 
This symmetry will be used extensively in the upcoming sections. 
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Proposition 2.1 Equation ( [ff. !\ ) remains invariant under 

(t,x,y,z) -> (-t,x,-y,-z). 

In order to understand the behavior described in the succeeding sections, it is first 
necessary to understand the dynamics of fl2.7|) when e = 0. In Figure [l] the flow on the 
invariant (when e = 0) plane {z = 0} is shown for different values of fj,. When /i < 
there exists a heteroclinic connection which connects the invariant plane {x = 0} to itself. 
This connection corresponds to a bright solitary wave. For /i = there exists a connection 
between x = and x = xq. This solutions corresponds to a kink. Due to the symmetry 



described in Proposition 2.1, there also exists a connection between xo and zero. Finally, 
when [i > there exists a solution connecting xq to itself, which corresponds to a dark 
solitary wave. Although it is not shown in the figure, this solution is actually embedded in 
the intersection of two two-dimensional manifolds in the phase space. 

In the upcoming sections it will occasionally be useful to have an analytic expression for 
the front (kink) which exists when e = /i = 0. This will be especially true when attempting 
to gain information on coefficients used in asymptotic expansions. Fortunately, due to 
Marcq et al [?] such an expression is available. 

Proposition 2.2 When e = /x = there exists a solution (x(t),y(t), z(t)) = (p(t),u(t), (/>(i)) 



to (2.7) which is given by 



Pit) = ^(l + tanh^i))^ 

u(t) = p'(t)/p(t) 
</>(t) = 0. 

In the above xq is the largest positive root of X(x) = and is given by 



2_ 3 
X °~ 4^ 



(recall that a ~ — 0.1). 



In the subsequent sections the notation p ■ t will be used to represent the trajectory of 
a point p = (x, y, z) under the flow generated by ( |2,7p . Under this convention, p ■ = p. In 
addition, given a point p = (x, y, z), if it is of interest to only study the behavior of one of 
the variables, say x, under the flow, the notation x ■ t will be used. 

When discussing the geometric objects in the following sections, the reader should con- 
sult Figure ^ to aid in the visualization. 

3 Geometric objects 

For this section assume that p = p(e), with fi(e) = 0(e 2 ). The validity of this assumption 
will be verified in Section 4. 

3.1 Flow near {x = 0} 

The goal of this subsection is to characterized the flow near the invariant plane {x = 0}. 
First, there exists a pair of critical points, say (0, y±(e), z±(e)), which satisfy the algebraic 
equations 

-y 2 + z 2 - A(0) + M (e) = 

-2yz - e(uj(x ) + ea) =0. { ' ' 
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It is not difficult to check that 

V±(0) = 

furthermore, since //(e) = 0(e 2 ), 

y±(e) 



2/±(0)+O(e 2 ) 



(3.2) 



It is also of interest to note that due to the symmetry described in Proposition 2.1, y_ = 
—y+ < 0, z_ = — z + . Furthermore, upon recalling the expression for xq given in Proposition 
2.2 , it is not difficult to see that A(0) = — Xq/4, and hence y±(0) = ±xo/2. 
For convenience, set 

7± = »±(0), <*± = -^y 
so that the critical points can be said to satisfy the asymptotic expansion 



y±(e) 

*fc(«0 



7± + 0(e 2 ) 
a±e + 0(e 2 ) 



Linearizing the vector field about the critical points (0, y±, z±) yields the matrix 

/ 



which has the eigenvalues 



A± 



4 




A4 








7± + 0(e 2 ) 

-2 7± -Ha ± e + 0(e 2 ) 



(3.3) 



(3.4) 



A 2 



Thus, the point (0, y+, 2+) has a one-dimensional unstable manifold, W"(0), coming out 
of the plane {x = 0} tangent to the vector (1,0,0), and a two-dimensional stable manifold 
embedded in the plane {x = 0}. In addition, the point (0,y_,z_) has a one-dimensional 
stable manifold, VF s (0), coming out of the invariant plane and tangent to (1,0,0), and a 
two-dimensional unstable manifold embedded in the invariant plane. 

Given 5 > 0, set 

C = {(x,y,z) : 0<x<5,y = z = 0}. (3.5) 

Due to the symmetry, this curve will be of paramount importance in subsequent sections. 
In the following proposition the term rj is assumed to be positive and small. 

Proposition 3.1 Let p = (x,y,z) G Co be such that x = 0{e n ) for some n > 1. Then as 
long as y_ + rj < y ■ t < y+ — 77, z ■ t 7^ for t ^ 0. Furthermore, z ■ (—t) = —z ■ t. 

Proof: By the nature of the function uj(x) it is clear that for x = 0(e n ), u(x) = 0(e 2n ). 
Since x' = xy, x ■ t = 0(x ■ 0) as long as y- + 77 < y ■ t < y + — -q. Thus, if x ■ = 0(e n ), then 
x ■ t = 0(e n ) for y ■ t in the prescribed range. By the above description of uj(x) for x small, 
this then implies that for x ■ = 0(e n ), 



-2yz - e(w(a? ) + ecr) + <3(e 
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If n > 1 it is now clear that trajectories can cross the plane {z = 0} at most once for 
y_ + 7] < y < y + — r/. This proves the first part of the proposition. The second part follows 
immediately from the symmetry inherent in Q2.7|). ■ 



It is now of interest to determine the nature of the curve Co as the flow carries it by 
the critical point (0, y+, z+). Using the matrix A + defined in ( |3.3| ), the linear flow near this 
critical point satisfies 



x' = ax 



y > = -2ay + 2ebz (3.6) 
z' = —2ebz — 2ay, 



where 



a = 7+ + 0(e 2 ), b = a + + 0(e). 
Setting y = rsinO, z = rcosQ equation ( |3.6| ) can be rewritten as 



x' = ax 

r' = -2ar (3.7) 

9' = eb, 

which has the solution 

x(t) = xoe at 

r(t) = r e- 2at (3.8) 

e\t) = Oo + ebt. 

For each p £ Co define to(p) to be such that 

t (p) = {inf : p ■ t n {y+ - r = r?} + 0}, (3.9) 

and set 

C -t = (J p-t (p). (3.10) 

After translating the critical point (0, y + , z + ) to the origin, the curve Co ■ to can be written 
parametrically as 

Co • to = {{x, r,9) : x = Xl (s), r = r,,9 = -vr/2 + ^(s)}, 

where xi(0) = 0, x[(s) > 0, xi(s) < rj, and #i(s) = O(e). 

It is of interest to determine the nature of Co • to as the flow forces it to intersect the 
plane {x = rf\. Given a point in Co • t , the time of flight, tf, is defined by x(tf) = r\. 
Sustitution of this expression into (|3.8|) yields 

tf = -ln^L-. (3.11) 

a xi(s) 

This expression for tf further yields that 

1 2, 



11 (3.12) 



2 a x\(s) 
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The above curve is a logarithmic spiral centered upon the point at which W M (0) first inter- 
sects {x = rf\. It is important to note that 

Xl (s) = 0(e n ) 9(t f ) = ~ + 9 1 (s) + 0(ne\n- e ) 

x 1 (s)=0(e-^) 0{t f ) = -- + e 1 {s) + O{l). 

Thus, the spiralling effect is seen only for those points which are initially exponentially close 
to the plane {x = 0}. 

Define the transverse sections to W u (0) and W s (0) by 



B o = ii x ,y, z ) x = r],\y -y^\ <r],\z - <7]} 

B o = {( x ,y, z ) ■ x = r i, \y -y+\ <v,\z - z+\ <7]}, 



where, as before, i] > is small. Now define 



= r(0)nB - = {(r],y-+0(r t ),Z-+0(r t ))} 
p u = W u (0)nB+ = {(7 1 ,y + +0( V ),z + + 0( V ))}, 



(3.14) 



(3.15) 



where the above represents the first intersection of the manifold with the set. 

Assume that 8 is sufficiently small so that for each p G Cq there exists a t > such that 
p-tnB^ ^9. For each p G C define 

tf(p)={inf : (3.16) 

and set tjT (p) = — i(j~(p). By using a symmetry argument it can be shown that p-t$ (p) G -B^. 
Finally, define the flow of Cq as it intersects the sets by 



C o • t = \J p-t (p) 

peCo 

Co -4 = U p-^)- 

pec 



(3.17) 



When e = these sets are given by 



Co-t = {{ x ,y, z ) '■ x = v,y- < y < y- + v,z = 0} 

Co -to = {{x,y,z) : x = rj,y + - rj < y < y + ,z = 0}. 



(3.18) 



The blow-up sets will be defined next. This set corresponds to points p such that 
\p ■ t\ — > oo in finite time ([?]). As a preliminary, set 



B 0,s = {( x ,y, z ) ■ < x < r/,y = y^ -r],\z - Z-\ <7]} 
B 0,s = {( x ,y, z ) ■ < x < r/,y = y + + r/,\z - z + \ <7]}. 

Following Kapitula and Maier-Paape [?] there exist sets 

C-oo = {P ^ B o~s '■ 3 t{p) > such that lim p ■ t = (0, -oo, 0)} 

t->t- (p) 

C+oo = {p£ #o~ s : 3 tip) < such that lim p ■ t = (0, +oo, 0)}. 

t->t+(p) 



(3.19) 



(3.20) 



The above sets are smooth curves, and symmetry considerations yield that (x, y, z) G C_ oo 
implies that (x, — y, — z) G C +oc . For p G C+oo set 

i +0O (p) = {inf : yinB+^l, (3.21) 
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and for p £ C_ ra define 

*-oo(p) = {sup : p-tnBo ^0}. (3.22) 

*<o 

Finally, the flow of these sets will be given by 



C-oo • t-oo = [J p-t-oo(p) 



(3.23) 



When e = these sets are given by 



-oo u — oo 



{(x, y, z) : x = rj, y_ - rj < y < y_, z = 0} 
C+oo-i+oo = : x = ??,y + < y < y+ - n,z = 0}. 



(3.24) 



This subsection can now be concluded with the following lemma, whose conclusion 
follows from the above discussion. The fact that C-oo • t-oo and C+oo • t +OQ are logarithmic 
spirals follows from the same argument as that which led to the description of Co • £q and 
Co- 1+. 

Lemma 3.2 Suppose that ew(xo) 7^ 0. Then Co • t^ and C-oo • t-oo are logarithmic spirals 
centered on Pq, and the curves Co ■ and C-oo ■ ^+oo are logarithmic spirals centered on 
Po- The spiralling is seen only 0(e _C//e ) close to the points p^ and p^; furthermore, their 



location is described within 0(e) by equations ( 3.18 ) and ( 3.24 )■ 



Remark 3.3 By the symmetry described in Proposition yl. i| it is necessarily true that 
(n,y,z) G Co • ifnplies that (rj, —y, —z) G Co • t$, and (n,y,z) G C-oo ■ t-oa implies 
that (r/, -y, -z) G C +oc • t +OD . 

3.2 Flow near M e 

The purpose of this subsection is to describe the flow generated by Q2.7| ) near the point 
(xo,0,0), where xo is the largest positive zero of A(x) = 0. 

When e = there exists a critical manifold A4o which is given by 

Mo = {(x, y,z) : y = 0, z 2 - A(x) = 0, Az 2 + xA'(x) < 0}. (3.25) 

This manifold is normally hyperbolic, and therefore smoothly perturbs to a manifold M t for 
e sufficiently small. Furthermore, the manifold M t has a two-dimensional stable manifold, 
W s (M e ), and a two-dimensional unstable manifold, W u (M e ), which are smooth perturba- 
tions of the center-stable and center-unstable manifolds which exist when e = (Fenichel 
[?], Jones [?]). 

As in the previous section, assume that \i = /i(e) = 0(e 2 ). The critical points for 
which lie on A4 e satisfy y = and 



A(x) - z 2 + 0(e 2 ) = 
oj(x) — u>(xq) — ea = 0. 



(3.26) 



A Taylor expansion about the point xo yields that 

1 

V(x ) 



x-x = e— ^— + 0(e 2 ). (3.27) 
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Substituting the expression in ( |3.27| ) into the first equation of ( |3.26| ) and taking a Taylor 
expansion yields 

X ' ' z 2 + 0(e 2 ) = 0, (3.28) 



oj'{x ) 

from which one gets that 



z2 = £ V(xo) 2 

U'(X Q ) 



This in turn yields 



2 = Y' (l0> ^ el/2 + O(e) ' (3 ' 30) 



Letting (x*(e), 0, ±z*(e)) G ,M e represent the critical points, it is now seen that these points 
have the asymptotic expansion 

x*(e) = x + ^ T e + O(e 2 ) 

3.31 

+ 0(6). 



Since A'(xo) < 0, in order for the above expressions to make sense the following assumption 
must hold. 

Proposition 3.4 If the parameters di,d,2, and a are chosen so that 

to'(xo)a < 0, 

i.e., 

(dx + 2x 2 Q d 2 )o- < 0, 



then there exist critical points on Ai t whose x and z components are given by ( 3.31 ) 



The next goal is to get an asymptotic description of A4 e within 0(e) of {z = 0}. A 
preliminary lemma is first needed. In the following, let the curve C XQ be given by 

C XQ = {(x,y,z) : y = z = 0, \x - x \ < u}. (3.32) 



Lemma 3.5 There exists an eo > such that if < e < eq, then C XQ n M t ^ 0. 

Proof: Clearly, A4 e intersects the plane {z = 0} nontrivially. Set 

{(x e ,y e ,0)} =M e n{z = 0}. 

In order to prove the lemma, it must be shown that y t = for e / 0. By the definition of 
Mo and the fact that M e is a smooth perturbation of M$ it is clear that 

lim y t = 0. 

In a sufficiently small neighborhood of A4 t the manifolds W s {M e ) and W u {M e ) satisfy 
the relation 

W s {M e )DW u (M € ) =M e . 
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This immediately yields that 

W s (M e ) n W u (M e ) n {z = 0} = {(x £ , y e , 0)}. (3.33) 
Suppose without loss of generality that 

lim (x e ,y e ,0)-t=(x*,0,±z*). 

t— >±oo 

By the symmetry of the ODE, if (x e , y e , 0) • t is a solution, then so is — y e , 0) • (— t). This 
immediately yields that 

lim (x € ,-y e ,0) -t = (x*,0,=fz*), 

t— >±oo 

so that 

{(x e , -y e , 0)} c w s (.M £ ) n W U (M £ ) n{z = 0}. 



But equation ( 3.33 ) then gives 

(x e ,y e ,0) = (x e ,-y e ,0) 

so that y e = 0. ■ 

For the rest of this subsection set ^ = \z\ + e. The slow manifold is given by the 
graph 

M e = {{x, y,z,e) : x = M x (z, e),y = M y (z, e)}. (3.34) 

By the definition of Aio the function Ai y satisfies the relation A4 y (z, 0) = 0. Furthermore, 
the conclusion of Lemma 3.5 yields that A4 y (0,e) = 0. Thus, a Taylor expansion of A4 y 
gives 

M y (z, e) = ezM y (z, e), (3.35) 

where the function A4 y is uniformly bounded for £ sufficiently small. The following lemma 
gives a description of the functions M x and My 

Lemma 3.6 The functions M x and M y comprising the graph of M e have the Taylor 
expansion 

M x (z,e) = x + -——z 2 + xe 2 + 0(C 3 ) 
\'{x ) 

where x is an O(l) constant. 

Proof: The basic idea of the proof is to write out a Taylor expansion for both M x and M y 
and to use the fact that the manifold M e is invariant under the flow. After using (3.35) 
the functions can be expanded as 

M x (z,e) = x + x\z + x 2 e + x 3 z 2 + x 4 ez + x 5 e 2 + 0(( 3 ) 
M y (z,e) = ez(y + yi z + y2e + O(( 2 )), 

where xq is such that X(xq) = 0. In order to prove the lemma, it must be shown that 

xi = x 2 = x 4 = y = y x = 
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and 

1 2ax 3 

When e = the function A4 X satisfies the relation 

z 2 - X(M x {z,0)) =0. 

Taking a Taylor expansion of A(x) about xq and equating coefficients quickly yields the 
desired relation for x\ and X3. 

In order to show that X2 = 0, consider the following argument. The flow on A4 e satisfies 

z' = -2M y (z, e)z + e(u{M x (z, e)) - u(x ) - &r), (3.36) 

so that when z = 

z' = e((x 2 co'{x ) - a)e + 0(e 2 )). 

The above expression is arrived at upon taking a Taylor expansion of tu{x) about xq and 
using the expansion for Ai x (0,e). Note that z' = 0(e 2 ). Now, the fact that on M e , y = 
M y (z,e) implies that 

y' = z'd z M y (z,e). 

Using the Taylor expansion for M y and the fact that z' = 0(e 2 ) yields that when z = 
0, y' = 0(e 3 ). But when z = 0, 

y > = -M 2 y (0, e) - \(M X (0, e)) + /i(e). 

Since /u(e) = 0(e 2 ), after taking a Taylor expansion of X(x) about xq one gets that 

y' = -x 2 \>(x )e + O(e 2 ). 

Thus, the fact that y' = 0(e 3 ) necessarily gives that x 2 = 0. 

Using the above results and equation ( |3.36| ) gives that the equation for z on M. e is 

z' = -ae 2 + (x3Uj'(x ) - 2y )ez 2 + x 4 uj' (x )e 2 z + x 5 u>'(x )e 3 + 0(C 4 ). 

Since x' = xy, on A4 e 

M' x = M x M y . (3.37) 

But 

M' x = z'd z M x 

= -2ax 3 e 2 z - ax 4 e 3 + 0(C 4 ) 

and 

M x M y = x y ez + xoy^z 2 + x y 2 e 2 z + 0(C 4 ), 
so that ( 3.37| ) becomes 

-2axse 2 z - crx 4 e 3 + 0(C 4 ) = x yoez + x yiez 2 + x y 2 e 2 z + 0(( ). 

Upon equating coefficients the desired result is reached. ■ 
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3.3 Flow near M e 

Recall that the flow on Ai e is given by ( 3,36| ). Using the result of Lemma 3J3 and Taylor 
expanding uj(x) around x$ then yields that on jW e the flow is given by 

z' = e(^^-ae + 0{^)) + 0{C% (3.38) 

Note that the above equation implies that for \z\ = 0(e), 

z' = -e 2 a + 0(e 3 ). (3.39) 

This estimate will turn out to be crucial for the estimates that follow. 

In order to determine the nature of the flow near A4 e , it will be desirous to use Fenichel 
coordinates ([?], [?], [?]). Before doing so, it will first be beneficial to understand the matrix 
arrived at upon linearizing about M e . Specifically, the manner in which the eigenvalues 
and eigenvectors vary as a function of z and e must be understood. 

Linearizing the flow (|2.7| ) about a point (x, y, z) £ Ai e yields the matrix 

V x 

A M , = 




-X'(x) -2y 2z j (3.40) 

l euj'(x) —2z 

which has the characteristic equation 

P( 7 , z, e) = ( 7 - y)((7 + 2y) 2 + 4z 2 ) + x(\'(x)(<y + 2y) - 2eJ{x)z) = 0. (3.41) 

The above equation has two solutions which are 0(1), say 7 ± (z, e), with 7 + = —7" and 
7- (0,0) = -V-zoA'(xo). 

In the following calculations the result of Lemma 3.£ will always be implicitly used. It 
is easy to check that 

P 7 ( 7 - 0,0) = -2x A'(x ) 

and 

p,( 7 - 0,0) =P e (0,0) =0, 

so that 

7 - (0,0) =77(0,0) = 0. 

The above follows immediately from the fact that 

j- = -P z /P 7 , j- = -P e /P 7 . 
Taking second derivatives gives 



P«(7~ 0,0) = -8V-x A'(s ) 
P M ( 7 -,0,0) = -2x o; / (xo). 

Since 

Tzz — Pzz / P~f ■> ^zt — Pze/Pyj 

this immediately gives that 

7^(0,0) 
7,1(0,0) 



yJ-X Q \'{XQ) 

uj'(xq) 



X'(x ) 

The proof of the following proposition is now complete. 
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Proposition 3.7 The O(l) eigenvalues of Aj^ e satisfy 

1) j + (z,e) =-j~(z,e) 

2) rM = -v-x„vw + ^=^= 2 - ^« + *• + 0( a, 

where 7 is an 0(1) constant. 

Although the following information is not necessary at the moment, it will be useful at 
a later time. The eigenvector v 7 - of Aj^i e associated with the eigenvalue 7" is given by 

/ - -2(j~ -y)z + exw'(x). T , . 

v 7 - = (^,7 -I/, — 7-5 • ( 3 - 42 ) 

7 + 2y 



Let v~- = (f 1, t>2, ^3) T - After using the expansions given in Lemma 3J3 and Proposition 3/7 
the following proposition is realized. 

Proposition 3.8 The components of the eigenvector v 7 _ are given by the expansions 

vi = x + - T }—z 2 + ie 2 + O(( 3 ) 
A'(xq) 

„ = -2,- ;"">;> e+0(c 2 ), 

where x and 7 are O(l) constants. 

Remark 3.9 Due to the symmetry of the ODE, the eigenvector associated with 7 + is such 
that the second and third components are the negative of those given above. 

Remark 3.10 When the above expressions are evaluated at the critical points (x*(e), 0, ±z*(e)), 
it is of interest to note that the third component, V3, satisfies the estimate 

v 3 = -2(±z*(e))+0(e). 

The above expression is valid because it is known that z*{e) = 0(e 1 ^ 2 ). 

Recall that for \z\ = 0(e) that the flow on A4 e satisfies 

z' = -eV + 0(e 3 ). 

Thus, in Fenichel coordinates the equations near Ai t for \z\ = 0(e) may be written as 

a' = A(a, b, z,e)a 

b' = T(a,b,z,e)b (3.43) 
z' = e 2 (—cr + h(a,b,z,e)ab). 

In the above equation the set a = refers to W s (Ai e ), the set b = refers to W u (A4 e ), 
and the function h is uniformly bounded. Using the eigenvalue expansion in Proposition 



3.7| one can say that for \z\ = 0(e), 



A(0, 0, z, e) = J-x X'(x ) + 0(e 2 ), T(0, 0, z, e) = -A(0, 0, z, e). 
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After taking a Taylor expansion the equations (|3.43|) may then be rewritten as 



a' = ^- x X'(x )a + O((\a\ + |ft| + e) 2 ) 

b' = - v /-x A'(x )6 + O((|a| + |6|+e) 2 ) (3.44) 
z' = e 2 (—a + h(a,b, z, e)ab). 

Define the set B XQ by 

B xo = {(a, b, z) : \a\ + \b\ < u, \z\ < O(e)}. (3.45) 

For e and v sufficiently small it is clear that one can, without loss of generality, use the 
linear equations 

a! = y/- x \'(x )a 

b' = -^- Xo X'(x )b (3.46) 
z' = -e 2 a 

when discussing the flow in B XQ . For convenience, in the rest of this subsection 7 = 
y/—xo\'(xo), so that ( |3.46| ) can be rewritten as 

a' = 7a 

b' = -76 (3.47) 
z' = -e 2 a. 



Now that the equations for the flow near Ai e have been determined, it will be of interest 
determine the nature of the set C Xo 
precise, for each p £ C XQ define t XQ by 



to determine the nature of the set C Xo as it exits B XQ under the influence of the flow. To be 



tx (p) = {sup : p-t € B+}, 

and set 

C X0 -t X0 = |J p-t X0 (p). (3.48) 
The set C XQ ■ t XQ will be close to the set a = 0, i.e., close to W s (Ai e ), as it exits B£ , where 



B t = {(x,y,z) : x = x - u,\y - (\J -\'(xq)/xq)v\ < v, \z\ < u} (3.49) 
The following lemma gives a determination as to how close. 

Lemma 3.11 Let p £ C XQ be such that t XQ (p) = 0(1/ e). The curve C xo -t xo is C 1 — 0(e~ c ^ e ) 
close to a = in a neighborhood of p ■ t xo (p), where c > is some constant. 

Proof: In order to prove the lemma the time-reversed flow for ( p. 47 ) must be considered. 
After such a time reversal, the solution to ( |3.47| ) is given by 

a(t) = aoe~ 7 * 

b(t) = 6 e 7 * (3.50) 
z(t) = z + e 2 at. 

In Fenichel coordinates the curve C XQ can be parametrically represented as (ao(s), bo(s), zo(s)) 
(s, s, 0) with < s < v. Substituting this set of initial conditions into ( 3.5C| ) yields that the 
flow of C Xo is governed by 

a(t) = se-^ 

b(t) = se^ (3.51) 
z(t) = e 2 at. 
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A solution to (3.51) leaves B xo when b(t) = v. Solving this equation yields that the 
time-of-flight, tf, is given by 



1 v 
t f = - In-. 

7 s 



Substituting tf into ( |3 . 5 1| ) yields 

«(*/) = 77 



,s 2 



z(tf) = -e - In-. 

7 s 

The above equation is a parametric representation of a curve on the section B XQ . It can 
be rewritten as 

a(z) = ue Pz/e \ 

where 

P = 2j/a, az < 0. 
The conclusion of the lemma is now clear, as when z = 0(e) with az < 0, 

a(z) = 0(e~ c A), a'(z)=0(^), - 



4 W u {0) n W s (.M e ) transversely 

After appending y! = and e' = to (^]) one arrives at a three-dimensional mani- 
fold W"(0) (abusing notation here) and a four-dimensional manifold W S (M 6 ) in the five- 
dimensional phase space. Since x' > along the relevant trajectory for < x < xq, by 
continuity this holds for e and \x sufficiently small. As such, for < x < x$ the manifolds 
can be parameterized in the following manner: 

W u (0) = {(x,y,z) : y = Uy(x,fi,e),z = U z (x,fi,e)} 

W s (M e ) = {(x,y,z) : y = Sy(x,z,v,e),\z\<e<l}. { ' ) 

A heteroclinic orbit possibly exists when W"(0) D W s (M t ) nontrivially, or equivalently 
when the function 

F(x, n, e) = S*(x, W(x, n, e),n, e) - W(x, fi, e) (4.2) 

has a zero. 

It is clear that F(x, 0, 0) = 0. If F^(x, 0, 0) / for some value of x G (0,xq), then 
an application of the implicit function theorem yields that there exists a smooth function 
fi = n(e) such that F(x, //(e), e) = 0. In other words, the condition that F^ ^ gives that 
the manifolds W u (0) and W^A^) intersect transversely in (x, y, z, //)-space, so that their 
intersection persists for small perturbation. Furthermore, the implicit function theorem 
yields that 

F £ (x,0,0) 

^ (0) = -i^o)- (43) 

Thus, it is possible to get asymptotics for the function //(e) by thoroughly understanding 
the manner in which the manifolds intersect. 
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Differentiating (4.2) yields that 



F^x, 0, 0) = (5* - Uj> + SJlU*)(x, 0, 0). (4.4) 

An examination of ( |2.7[) yields that when e = it is unchanged under the transformation 
z — > —z; therefore, 

5*(z,0,/i,0)=0. (4.5) 
This gives the simplified expression of F^ to be 

F^x, 0, 0) = S*(<c, 0, 0, 0) - W*(s, 0, 0). (4.6) 

In a similiar manner, it is seen that 

F e (x, 0, 0) = S y e (x, 0, 0, 0) - U%(x, 0, 0). (4.7) 

Let the heteroclinic orbit which exists for p, = e = be denoted by (p(t), u(t),0). Assume 
that the wave has been translated so that p(0) = xi, where < x\ < xq is the other root of 
X(x) = 0. When linearizing ( |2.7D about this orbit the variational equations take the form 

Sx' = u5x + p5y 

by' = -\'{p)5x - 2u5y + 5p 

Sz' = -2u5z + (u(p) - u(x ))8e (4.8) 
Sp! = 
5e' = 0. 

Each of the coordinates 5xi of ( |4.8| ) (x{ = x, y, etc.) is a 1-form. From these 1-forms exterior 
products of forms of any degree k can be constructed. The rest of the following discussion 
will only be concerned with 2-forms P XiX = 5xi A5xj; however, it can be generalized to any 
k. 

Each 2-form associates to a 2-plane T a number that is the area of the projection of 
a unit square of T onto the coordinate planes of the two coordinates specified by Px^Xj ■ 
In order to evaluate P XiX j, the following rule is applied. Let ./V represent a /c-dimensional 
manifold, k > 2, and let T p N represent the tangent space to N at a point p G N. Let 
{ a i(p)} represent a basis for T p N. Suppose that S(p) C T p (n) represents a two-dimensional 
subspace which is spanned by the vectors {a,i(p), a,j(p)} for some 1 < i,j < k. These vectors 
may be thought of as rows to a 2 x k matrix, A(p). If Xi and Xj represent two of the k 
coordinates of T p N, then P XiX - is the determinant of the 2x2 submatrix obtained by looking 

at the x^ 1 and columns of A(p). Furthermore, the evolution equation for P XiXj is given 
by the product rule, i.e., 

P' XiXj = 5x' i A 5xj + 5xi A Sxj. 
Lemma 4.1 For any < x < xq, U^(x,0,0) > 0. 

Proof: Set P xy = Sx A Sy and P Xfl = Sx A 5p. Using the equations of variation ( |4.8j ) it is 
seen that 

^xy — uP X y + uP x ^. 

Two vectors which are tangent to the manifold iy u (0) at a fixed x value are 

ei = (//y,o,o,o) T , 6 = (o,^,o, i,o) T . 
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In the above the fact that U z (x, p,, 0) = is implicitly used. When applied to these vectors 

so that the equation for P xy becomes 

p 'xy = -uP xy + pu 2 . 
Since u = p' / p, this equation has the solution 

p{t)P xy (t) = [ p 2 (s)u 2 {s)ds, 

J —oo 

so that P X y{t) > for all t. Since for a fixed x value 

P^(£i,£ 2 )ocp'W#M,0), 
the fact that p' > yields the conclusion of the lemma. ■ 

Lemma 4.2 For any < x < x , 0,0,0) < 0. 

Proof: The proof follows that of the previous lemma, and will therefore only be sketched. 
Set 

6 = (p>',0,0,0) T , 6 = (0,5^,0, 1,0) T . 
The equation for P xy , when applied to these two vectors, becomes 

Ky = ~uP xy + pu 2 , 

which has the solution 

xy 

Since for any fixed x 

Pay (6, 6) 0Cp'5»(x, 0,0,0), 
the conclusion immediately follows. ■ 

Corollary 4.3 For any < x < xq, F fl (x,0,0) < 0. 

Proof: The result follows immediately from the above lemmas and equation ( |4.6| ). 
Lemma 4.4 U]/(x,0,0) = 0. 



oo 

p(t)P xy (t) = - I p 2 (s)u 2 (s)ds. 



Proof: The proof follows that of Lemma 4.1, and will therefore again only be sketched. 
Set 

a = GoV,o,o,o) T , 6 = (0,^,0,0, i) T . 

When applied to these two vectors, the equation for P xy is 

Pxy = ~ u Pxy, 
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which has the solution 

Pxy{t) = WY 

Since limt^-oo P xy {t) must be bounded, the fact that p(t) — > as t — ► — oo implies that the 
constant C must be zero. Since for fixed x 

P*»(6,6)oc^(x,0 J 0) J 

the conclusion of the lemma now follows. ■ 

The proof of the next lemma mimics those of the previous three lemmas, and will 
therefore be left to the reader. All that must be kept in mind is that the quantity P X y(t) 
must approach zero as t — ► oo, as £i(0 — ► and the vector £2 CO remains bounded as t — ► 00. 

Lemma 4.5 5* (a, 0,0,0) = 0. 



Corollary 4.6 F £ (x,0,0) = 0. 

Proof: The conclusion follows immediately from the previous two lemmas and equation 



As a consequence of Corollaries |4.3| and |4.q , this section can be closed with the following 
theorem. 

Theorem 4.7 There exists a smooth function // = //(e) such that when < e <C 1 i/ie 
manifolds W s (M. e ) and W u (0) have a nontrivial intersection. Furthermore, the function 
//(e) satisfies 

//(0) = M '(0) = 0, 

so f/iaf 

//(e) = 0(e 2 ) 



Remark 4.8 Due to the symmetries present in (2.7), when // = //(e) it is necessarily true 
that W s (0) Pi W u {M t ) is also nontrivial. 

5 W u {M e ) n VK s (>( e ) transversely 
5.1 Flow from M. e to {x = 77} 

Assume that // = //(e) = //2e 2 + 0(e 3 ). Prom the results of Section 4 it is then known 
that W u (0) n W s (M e ) is nontrivial, as is W s (0) D W u (M t ). Furthermore, on M e there is 
a connection between the critical points (x*(e), 0, ±z*(e)). Since there exists a connection 
between (0, y_, z_) and (0, y + , z + ) in the plane {x = 0}, it is clear that when // = //(e) there 
is a heteroclinic cycle in the three-dimensional phase space. The goal of this section is to 
explore the possible solution set as // is varied from //(e). 

Initially, the focus will be on the intersection of iy u (0) and W s {M e ); however, this is 
not a restriction in any way, as the symmetry described in Proposition |2.1| immediately 
yields similiar results for the intersection of VF s (0) and W u {M e ). The heteroclinic wave is 
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singular; however, it can be written in a regular perturbation expansion for < x < xq — v. 
To set notation, for < x < xq — v write the wave as 

pit) = Po (t) + O(e) 

u(t) = u (t)+O(e) (5.1) 
4>{t) = eMt) + 0(e 2 ). 

The following lemma gives a characterization for <p\. 
Lemma 5.1 For t E (—00,00) the term (pi is given by 

Mt) = Pi{t)di + (3 2 {t)d 2 , 

where 

(3i(t) = -x 

(3 2 (t) = —4(5 + tanh(^i)). 

Remark 5.2 Note that 

(3i(t) 1 
f3 2 (t) 2x1 

for each t. 



Remark 5.3 Although it is permissable to use <t>\ when t G (—00, T v \, where p(T v ) = xq — v, 
note that f3\(t) and (3 2 (t) are well-defined for all t. 



Proof: After substituting the expansion ( |5.1[ ) into the ODE (2.7) and equating terms it is 
seen that <f>\ satisfies 

(f/ x = -2u <j)i + u(po) - co(x ). 
Since no = p'q/pq, the above equation clearly has the solution 

Po(. t )Mt)= I pl(s)(u(po(s)) - u(x )) ds. 
J — 00 

Thus, after using that to(a) = d\o? + d 2 a one finds that 

Mt) = Pi(t)di + (3 2 {t)d 2 , 

where 

m) = ^ ) f_jl{sM(s)-xi i )ds. 



Substituting the expression for po(t) given in Proposition 2.2 into the above expression gives 
the result of the lemma. ■ 

Corollary 5.4 When e 7^ the function 4>(t) satisfies 

$(T V ) = -x (d 1 + x 2 d 2 + 0(n))e + 0(e 2 ) 
0(T„) = -xo^ + ^xfa + O^e + Oie 2 ). 
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Proof: The result follows immediately from the expansion of the wave given in (5.1) and 
the above lemma. ■ 

It is of interest to understand the behavior of W s {M e ) as it intersects the section Bq. 
For fixed e the intersection of W s (M e ) with Bq forms a curve y = S y (r],z,e). Setting T n 
to be such that p(T v ) = 77, it will be desirable to determine 5f (77, e). In order to find 

this quantity, the tangent space to the manifold must be tracked. 

In the four-dimensional phase space the manifold W s (A4 t ) is three-dimensional and can 
be written as the graph 

W s (M e ) = {(x, y,z,e) : 77 < x < x , y = S y (x, z, e)}. (5.2) 

Over the underlying wave, for each fixed x the manifold's tangent space is spanned by the 
three vectors 

6 = (p>',0',O) T 

6 = (0,Sf,l,0) r (5.3) 
£3 = (0,Sf,0,l) T . 

Since W s {M e ) is a three-dimensional manifold, it will be necessary to use three-forms to 
properly track its tangent space. In particular, the forms 

Pxye = fix A 5y A fie 

and 

Pxze = fix A 5z A 5e 

will be used. Note that for fixed x 
and in particular 

5^^(T,),6)oc%^. (5.4) 
P \- L v) 

Thus, it will be desirable to compute P xye at x = 77. 

After linearizing about the wave for e / the variational equations become 



5x' = ufix + p5y 

6y> = -\'(p)5x - 2u5y + 2<j)8z + (2/j 2 e + 0(e 2 ))5e 

5z' = euj'(p)5x — 2cj)5y — 2u5z + [to{p) — oj{xq) — ea)5e 

fie' = 0. 



(5.5) 



The ODE for P xye is then 

P' xye = -uP xye + 2(pP xze , (5.6) 
which, using the fact that u = p' / p, has the solution 

p{t)P xye {t) = (x - u)P xyt {T v ) -2^ p(s)<P(s)P xze (s)ds. (5.7) 

It is clear that the evolution of P xze must be understood in order to understand the 
evolution of P xye . Towards this end is the following proposition. 



22 



Proposition 5.5 When e = 0, 



2 «o(t) 



Proof: In order to prove the proposition, it will first be shown that the vector £2 is given 

by 

&(t) = (0,0,x 2 /p 2 (t),0f. 

Recall that the vector £2 is formed by taking the derivative of the manifold W s (M e ) with 
respect to z. Since the flow has the symmetry (x, y, z, t) — ► (x, y, —z, t) when e = 0, it is 
clear that the first two components of £2 must be zero. From the variational equations the 
third component satisfies the ODE 



Sz = —2uq5z, 



which has the solution 



5z(t) = C/p 2 (t). 
The constant C can be chosen so that 



lim 6z(t) 



t— >oo 



1, 



from which arises the characterization of fe(t)- 

Since = (p' (t),u' Q (t),0,0), the conclusion of the proposition now immediately fol- 
lows, as 

* 

r-2 /„2 



xl/ P l * 

1 



Recall that T v and T v are defined by 

p(T v ) = rj, p{T u ) = xq v. 
For t £ [T V ,T U ] the three-form P xze is given by the regular perturbation expansion 



2 u (t) 



P xze (t)=x z ^y- + O(e). (5.8) 



Substitution of this expression into ( |5.7|) along with the expansion for (f> given in Lemma 
5.1 yields that 

Po(t)P xye (t) = {x Q -v)P xye {T u )-2exl J^\ {s){(5 1 {s)d l +I3 2 {s)d 2 )ds + O{e 2 ). (5.9) 

In order to finish the calculation, P X y e (T v ) must be determined. 
As it has already been seen, 

P xy e{T v ) = p'{T u )Sy(x Q - is, <j>{T v ), e); 

thus, the quantity S^{xq — v,(f)(T v ),e) must be calculated. Since v is small, this can be 
determined by understanding the variation of the linear approximation to the manifold 
W s (M e ). Since M e C W s (M e ) it is necessarily true that 

Sy(M x (z,e),z,e) =M y (z,e). (5.10) 
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Recall that Lemma 3.6 states that the manifold A4 e satisfies 

d z M y (z,e) = 0(( 2 ). 

As such, due to ( |5.10D the slow manifold has no effect on W s (Ai e ), at least up to 0(e). 
Therefore, in order to understand the variation of W s (M t ) it it enough to determine how 
the vectors which are tangent to W s (A4 e ) at A4 e vary with respect to z and e. Recall 
Proposition 3.5, which states that the tangent vectors to W s (M e ) at A4 e are given by 



where the quantities V{ are defined in Proposition 3.8. In order to determine SJf, it is then 
enough to calculate d z (v 2 /vi), where these quantities are defined in Proposition 3JS. A 
simple calculation yields 



v 2 uj (x ) 2 
-dz— = e— — — + 0(C ), 

Vl XqX'(Xq) 



from which it is seen that 



SKxo - v, #T„), e) = + OH) + 0(C 2 )- 

Since p'(T u ) = 0(u) one then has that 

Pxye{T u ) = 0(v)e + 0(C 2 ) 

Upon combining the above results, and using the fact that Pi(t) is bounded for i = 1,2 and 
for all t, it is seen that 

Po{t)P xye {t) = (ai(i)di + a 2 {t)d 2 + 0{u))e + 0(e 2 ), (5.11) 

where 

roc 

ai(t) = -2xlj u (s)(3 i (s)ds (5.12) 

for i = 1,2. 

Proposition 5.6 When evaluated at t = T v , 

ai(T„) = —2x0 In — 

Xo 1 2 
a 2 (T v ) = -2x 5 ln^ + -x 5 (l-^). 

Xq Z Xq 

Proof: The expression for ai(T„) only will be proved, as the proof for a 2 (T ri ) is similiar. 
For the rest of this proof set k = xq/2. 

Using the expression for po(t) given in Proposition ||]2] it is not difficult to see that 

Mt) = Po(t)/Po(t) 

xq sech 2 (kt) 



4 l + tanh(H) 
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Since, by Lemma 5.1, 0i(t) = — xq, it is clear that 



. / sech 2 (s) 



kt 1 + tanh 



= 4(in(2) - ln(^). 

Xq 

The last part of the above expression comes from the fact that 

1 + tanh(fct) = 2pl{t)/xl. 
Since po(T n ) = rj, upon simplification the expression for a\{T v ) follows. ■ 
Remark 5.7 Note that 

„ (rp \ > X Ql 

with 

ai(T v ) 4 x 

for < -q <C 1 . 

Now that ai(T„) has been calculated for i = 1,2, the following lemma has essentially 
been proved. All that is now necessary is to recall (|5.4|), the fact that < p'(T v ) = 0(rj), 
and the fact that y + = xq/2. 

Lemma 5.8 The manifold W s (M e ) satisfies the estimate 

S%(r], 4>{T V ), e) cx —z(aid\ + a 2 d 2 + 0(y))e + 0(e 
if 

where 

a\ = — Axq In — 

X ° 2 

a 2 = -Axnln — + Xnfl - -9). 

x xi 



2\ 



Remark 5.9 Note that a\/a 2 satisfies the same estimate as that given in Remark 5.7. 



Corollary 5.10 The manifold W u {M e ) satisfies the estimate 

W z (r J -<t>(T v ),e)=$y(r),<t>(T v ),e). 



Proof: The result follows immediately from the symmetry outlined in Proposition 2.1. 
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5.2 Flow from {x = rf\ to {y = ?/_ + i]} 

Recall the transversality argument given in Section 4. It was shown how the manifolds 
W u (0) and W s {M e ) vary with respect to // along the wave when e = 0. Due to the 
smoothness of the flow, these estimates persist for e and // sufficiently small. Furthermore, 
as a by-product of the symmetry discussed in Proposition |2.1| , the manner in which VF s (0) 
and W u (Ai e ) vary with respect to // is also understood. 
Specifically, the following is known. Set 

W u (M e ) n {x = V } = {(y t z,ii,e) : y = U»(z, //,e)} 

W s {0)n{x = r,} = {(y,z,//,e) : y = ^(//,e),z = £*(//, e)}. { °' ° } 



As a consequence of Proposition |2.1| and Lemmas |4.1| and [4.2| it can be concluded that 

5«(0,0) <0, W» (0,0,0) >0; (5.14) 

furthermore, these estimates hold for e and // sufficiently small. When // = //(e) the mani- 
folds intersect, so that 

5"( M ,e) = Z^(S z (//,e),//,e). (5.15) 

Thus, when e is fixed, by the estimate fl5.14[ ) for // > //(e) the manifold W s (0) is below 
W u (Ai e ), while for // < //(e) the configuration is reversed. 

Recall Lemma ^2, which gives a description of the curves Co ■ and C±oo ■ t±oo- Specif- 
ically, it is of interest here to contemplate the nature of Co • 4q and C_oo • t-oo both of which 
are contained in Bq and centered upon 

P s = W s (0)nB . (5.16) 

The set 

C = C-oo • t-oo U Pq U C ■ to (5-17) 

defineds a curve in Bq-, in fact, for eu;(a:o) 7^ it is a two-armed logarithmic spiral centered 
upon^Q. When fi = /u(e), W u (Ai e ) intersects Cq at p^; indeed, the intersection is transverse. 
Thus, as \x is varied from //(e), W u (M e ) continues to intersect Cq] however, the intersection 
no longer occurs at p Q . This implies that W U (M £ ) intersects either C_oo • t_oo or Co • tg~ . 
From ( p. 14 ) and Lemma 3.2 it can be concluded that if |/z — //(e) | is 0(e n ) for n > 3 but is 



not exponentially small in e, then 

//>//(e) => W u (^ e )nC -t #0 , , 

//<//(e) ^ W«(7W e ) n C_oo • t_oo 56 { °- L0) 

(see Figure |H). Furthermore, there exists only one point in each of the intersections. If 
< I//— //(e) I < 0(e -c / e ), then there exists a finite number of points in both W u (A4 e )riCo-t ~ 
and W u (Ai t ) n C_oo • t_oo, with the number of points increasing to infinity as // — ► //(e). 

For the rest of this subsection it will be assumed that |/z — //(e) | is not exponentially 
small in e. This is done to clarify the following arguments. In addition, it will be assumed 
that // > //(e); however, this is not a necessary restriction and the below arguments can be 
modified to draw conclusions if // < //(e). 

Since // > //(e), by ( 5.18| ) W u (M. e ) D Co ■ tp is nonempty, which immediately implies 



that W u (A4 e ) n Co is nonempty. Due to symmetry outlined in Proposition 2.1 this yields 
not only that W s (M. t ) fl Co is nonempty, but 

W u (M e ) n C = W s {M e ) n C C VF"(7W e ) n W s (M e ). 
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Thus, there is a connection between (x*(e), 0, z*(e)) and (x*(e), 0, — z*(e)) which passes 
"near" (to be made more precise later) the plane {x = 0}. The goal of the rest of this 
subsection is to understand the nature in which W u (M e ) and W s (Ai e ) intersect at Cq. 

Now that the orientation of the manifolds W u (M e ) and W s (Ai e ) is understood as they 
intersect the plane {x = rj}, it is desirable to understand their behavior under the flow as 
they pass near the critical points (0, and (0, y+, z + ), respectively. The goal of this 

subsection is to show that under a suitable restriction the passage near these critical points 
does not effect their orientation. It will be sufficient to track the manifold W u (M e ), as the 



symmetry present in (2/7) allows one to then draw an immediate conclusion regarding the 
behavior of W S (M € ). 

First, suppose that fj, is such that 

0<O(e" +1 ) < \fi-fi{e)\ <0{e n ) (5.19) 

for some n > 3. As it will be seen, this restriction on /i will guarantee that the spiralling 
action near the critical points on {x = 0} will be unseen by W u (A4 e ) as it passes near 
{x = 0}. Since n > 2, the perturbation of \i given in ( |5.19| ) will not effect the asymptotic 
estimates given in Lemma 5^ and Corollary [5.10 . Using the asymptotics for U y , it will be 



desirous to look at the image of the curve 

x = r], y = y + ecz, (5.20) 

where \z\ <r],yo = 0(e n ), and 

sgn (c) = sgn (U%). 

The point yo describes the difference between \i and //(e), and by supposition is therefore 
positive. 

It will be advantageous to look at the flow near the critical point in polar coordinates. 
Upon setting 

y = r sin 9, z = r cos 6 
the curve ( |5,20| ) can be parametrically described as 



r (s) = \]yl + lets + (1 + e 2 c 2 )s 2 

a n t -l/^o + ecs (5-21) 

e {s) = Tan ( ), 

s 



where |s| < rj. The linear flow near (0, y_,z_) satisfies the ODE 



x' = —ax 



r' = 2ar (5.22) 
& = eb, 



where 



a = ^=A(0) + O(e 2 ) 



Thus, upon using ( |5.21[) as an initial condition to (|5.22 ) it is seen that the behavior of Wf 
is well approximated by the solution formula 

x(t) = r]e~ at 

r(t) = r (s)e 2at (5.23) 
0(t) = d (s) + ebt. 
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It is now of interest to determine the nature of (5.21) as the flow forces it to intersect 
the cylinder {r = rf\. The time of flight, tf, is defined by r(tf) = r\. Using ( |5.23| ), this 
yields that 

-lln. V 



t 



f 



2a r (s) ' 

Thus, under the linear flow ( |5.23 ) there is the mapping from {x = rf\ to {r = rf\ given by 




y/Vy/ r o(s) 



V 



)(«) + e 7T ln 



V 



(5.24) 



2a ro (s) / 

for |s| < rj. Set x(s) to be the first component and 9(s) to be the third component of the 
vector on the right hand side of the above equation. Given the asymptotic expansion for 
it is then of interest to perform an expansion of x(s) and 9(s) around s = 0. This step 
is taken because it is only of interest to determine the behavior of the curve near where 
W u (A4 e ) intersects Co • t$. First, a Taylor expansion yields that 

ro(s) = yo + ecs + 0(w 2 ) 

6 (s) = l-± s + 0(w 2 ), 
2 yo 



where w = es/yo and it has been used that yo > 0. Plugging this expansion into (|5.24[) and 
performing another Taylor expansion then gives that 

ec 



V ' 2 2a y y K 2b 



2aC ^ s + 0(w 2 ). 



In rectangular coordinates the right hand side of ( |5,24 ) is given by x(s) 
and z(s) = rj cos 0(s). Using the expansion ( |5.25| ) it is then seen that 



y( s ) 



(5.25) 
77 sin 6>(s), 



x(0) 


= Vvy/m, 


x'(0) 


y(o) 


= r;sin^(0), 


y'(o) 


z(0) 


= rj cos 0(0), 


z'(0) 



ec 



2Vm 



r/(— (l + e 2 — )) cos 0(0) 



yo 



2b' 



(5.26) 



7? (I(l +e 2^ ))sin0(O) . 
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It is clear that the curve given by the right hand side of (5.24) can be parameterized by 
z for s sufficiently near zero. Using the expansion coefficients given in fl5.26| ) it is then seen 
that x = x(z) satisfies the relation 



dx 



dz 2^(l + e 2 §f))sin0(O)' 



(5.27) 



Since 



„(0) = f + £ Aln' 
2 2a y 



this means that for yo = 0(e n ), 



9(0) 



— h 0(ne In — ' 
2 v e' 



Substituting this estimate into ( ^.27| ) and applying ( |5.26| ) yields the following lemma. 
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Lemma 5.11 Suppose that fx — /x(e) = 0(e n ) for some n > 2. The intersection ofW u (M. e ) 
with the plane {y = y_ + rf\ is a parametric curve given by (x, z) = (x(s), z(s)). The curve 
satisfies the estimates 



b. 



c. 



x(0) 
dx 

dz 



0{e n ' 2 ) 

= 0( e (™+ 2 )/ 2 ) 



dx 

s 9 n (fa) = s 9n(Uy(r],-<j)(T v ),fi(e),e)). 



Remark 5.12 If W M (Ai e ) intersects the plane {y = y_ — rf\, then the conclusion of the 
lemma holds, with part c. being changed to 



sgn( 



dx, 
dz' 



-sgn(UV{r}, -<f>{T v ), /x(e), e)). 



As an immediate consequence of the symmetry of the ODE one gets the following corol- 
lary 

Corollary 5.13 Suppose that fi — ji{e) = 0{e n ) for some n > 2. The intersection of 
W S (M € ) with the plane {y = y + — rf\ is a parametric curve given by (x,z) = (x(s), z(s)). 
The curve satisfies the estimates 



a. x(0) = 0(e n / 2 ) 

b. — = 0( e («+ 2 )/ 2 ) 
dz s=0 

dx 

c- S9n<y -~aJ^ = s ^(5|(r/,0(T, ) ),/i(e),e)). 



Remark 5.14 If W u (M. e ) intersects the plane {y = y + + rf\, then the conclusion of the 
lemma holds, with part c. being changed to 

dx 

5.3 Flow from {y — y_ + i]} to {y = 0} 

Now that the passage of the manifold W u (Ai t ) near the critical point (0,y_,z_) is un- 
derstood, it is necessary to understand its behavior under the flow near {x = 0}. The 
previous lemma states that the manifold is 0(e n / 2 ) near {x = 0} when it intersects the 
plane {y = y_ + rj}. Since x' = xy and y < in the region of interest, it is clear that 
manifold stays within 0(e"/ 2 ) of {x = 0} until it intersects the plane {y = 0}. Thus, it is 
expected that the flow on {x = 0} dominates the behavior of the manifold. The goal of this 
subsection is to make this intuition rigorous. 

After the manifold W u (M e ) has intersected the plane {y = y_ + rj} it can be written 
as the graph 

W u (M t ) = {(x, y,z) : x= U x (y, z, e), y_ + r? < y < y+ - f], \z\ < rj}. (5.28) 

The previous lemma gives an expression for along the wave when y = y_ + rj. The goal 
is to calculate Wf(0, 0, e). In order to accomplish this, it will be necessary to once again use 
differential forms. 



29 



The vector to be tracked has the initial condition 

e= (U*(y- + r ] ,z,e),0,l) T , 

where the z-coordinate is taken to be such that the vector is over the underlying solution. 
The variational equations are 

Sx' 



Sy' 
5z' 

Projectivize the equations by setting 



uSx + pSy 
-\'{p)5x - 
euj'{p)5x 



- 2u5y + 24>5z 
2<p5y — 2u5z. 



5x/5z, b = 5y/5z. 



(5.29) 



(5.30) 



It is not difficult to check that 5z(£) = 0(1) for y_ +77 < y < 0, so that the above quantities 
are well-defined in the region of interest. The equations for a and b are given by 

a' = (—it + 2(pb)a + pb — euj'(p)a 2 
b' = -X'(p)a + 2(j)(l + b 2 )-eLo'(p)ab, 

with the initial conditions being 



(5.31) 



a(0) = U%(u-Q,<f)-0,e) 
6(0) = 0. 

Without loss of generality, assume that the z-coordinate of ^(Ojflfl^, —<j}(T v ), is such 



that — </>(T„) > 0. By Corollary 5.4, this implies that 

d\ + xld 2 > 0. 



(5.32) 



Let (p,u,(f)) -t £ W u (M e ) n W s (M e ) represent the orbit connecting (x*(e), 0, ±z*(e)), and 
let it be normalized so that u ■ = y_ + r/. It is known that this orbit intersects Co, as 
W u (A4 e ) n Co • i s nontrivial. Therefore, there is a To > such that u ■ Tq = (p ■ Tq = 0. 
In addition, there exists a T\ < such that p-T\ = 77. 

The claim is that 4>-t > for t £ [Ti, To]. To see that this is true, consider the following 
argument. The z-coordinate of W u (M. e ) n Co • £q > ' i s within 0(e 3 ), n > 3, of — <^(T^); 
therefore, since —^(T^) > is O(e), (f> ■ T\ is also positive. Following the argument of the 
previous subsection, it is not difficult to see that (f> ■ t > for t G [Ti, 0]. Since </> • > 0, 
as a consequence of Proposition 3.1 it is necessarily true that (j) ■ t > for t G [0, To] (recall 
that, by Lemma |ll|, p • = 0(e™/ 2 )). 

As a consequence of Corollary 5.10 and Lemma pMl it is known that 

sgn (Uz(u ■ 0, 4> ■ 0, e)) = sgn (aidi + a 2 c?2) 
W|(«-O,0-O,e) = 0( e ( n + 2 )/ 2 ). 

In what follows, assume that 

a\d\ + 0:2^2 > 0. (5.33) 

Since 01/02 is close to 1/xq (for an exact description see Lemma |5.8| ), this is not that much 
more a restriction than that already imposed by (5.32). 

Under this hypothesis, it is then being assumed that a(0) > 0, with a(0) = 0(e^ n+2 ^ 2 ). 
It is of interest to calculate a (To), as 



Wf(0,0,e) oc a(T ) 
Towards this end is the following lemma. 



(5.34) 
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Lemma 5.15 a(t) < 0(e n l 2 ) for t E [0,T ]. 
Proof: First, since p' = pu, 

p(t)=p(Q)efo u & ds , 

so that p(t) = O0(0)) = 0(e n / 2 ) for t E [0,T ]. In addition, since 8z(£) is O(l) for 
y_ + ?7 < y < 0, and <5x(^) and <5y(£) satisfy the linear ODE ( 5.29 ), it is necessarily true 
that a and b are 0(1) for t E [0,To]. These two facts will be used extensively in what 
follows. 
Set 

f(t) = e Sl u ( s )- 2< t> < y s ) b ( s ) ds 

_ P{t) -2 [*<P(s)b(s)ds 
P(0) 

Since b(t) = 0(1) and p{t) = O(p(0)), it is clear that f(t) = 0(1) for i E [0,T ]. Solving 
(|5.3l| ), it is seen that a(t) is given by the integral formula 

f(t)a(t) = o(0) + / f(s)p(s)b(s) ds-e t f {s)uj' {p(s))a 2 {s) ds. (5.35) 



Using the estimates on / and p, and assuming that a and 6 are O(l), it can then be 
concluded that 

f(t)a(t) = a(0) + 0(e n/2 ) • i + 0(e (n+2)/2 ) • t. 
Since a(0) = 0(e( ra+2 )/ 2 ), the conclusion follows. ■ 

It will be desirable to revise the above estimate. In order to do so, the following lemma 
is needed. 

Lemma 5.16 b(t) = 0(e) for t E [0,T ]. 

Proof: Since 6(0) = 0, the function b(t) satisfies the integral equation 

b(t) = 2 [ 0(a)(l + b 2 {s)) ds - [ \\' \p(s))a(s) + euj' {p(s))a(s)b(s)) ds. (5.36) 



Assuming that b(t) = O(l), which was justified in the proof of the above lemma, using the 
conclusion of the above lemma, and using the fact that <p(t) = O(e), one gets that 

bit) = O(e) -t + 0(e n ) -t. 

Since To = 0(1), the conclusion now follows. ■ 

A more careful examination of (5.36) yields the following result for b(t), which in turn 
can be used to improve upon the estimate made upon a(t) in Lemma 5.15| . Using the 
estimates provided in the above two lemmas, it is clear that 

b' = 2<p(l + b 2 ) + 0(e n ), 

which can be solved to get 

b{t) = tan( / 20(a) ds) + 0(e n ) • t. 
Jo 
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Substituting this expression into ( |5.35| ), and again using the above two lemmas, yields 
f(t)a(t) = o(0) + f f(s)p(s) tan( f 2<j)(r) dr) ds + 0(e 3n / 2 ) • t. 



Note that the second term on the right-hand side of the above equation is 0(e( n+2 )/ 2 ), as 
is a(0). Furthermore, this term is positive, as <f>(t) > for t £ [0, To]. Since, by supposition, 
a(0) > 0, it is now clear that 

f(T )a(T ) > 0(0), 

and furthermore, a(To) = O(a(0)). 

Using ( 5.34| ), the proof of the following lemma is now complete. While the lemma was 
proved only for the case of both <j> • > and a(0) > 0, it can easily be modified in the 
event that both quantities are negative. 

Lemma 5.17 Suppose that < \i — //(e) = 0(e n ) for some n > 3. Further suppose that 

sgn(ct\di + 02^2) = sgn(d\ +£0^2), 
where the on 's are defined in Lemma 15. <5|. Then 



sgn(U% (0, 0, e)) = sgn{a\d\ + a^cfe) 
Wf(0,0,e) = 0( e ( n+2 )/ 2 ). 



By the symmetries present in (|2.7[) one arrives at the following corollary. 

Corollary 5.18 Let the hypotheses of Lemma \5. 1 \ be satisfied. Then 

sgn(S*(0,0,e)) = -sgn(U*(0, 0, e)). 

Using the time of flight estimates given in Subsection 5.2, one finally gets the following 
theorem. 

Theorem 5.19 Suppose that < fi — /i(e) = 0(e n ) for some n > 3. Further suppose that 
the parameters d\ and d<i are chosen so that 

sgn(a\di +0^2) = sgn{d\ +xld 2 ), 



or where the ati's are defined in Lemma \5. <jj . Then the manifolds W u (M. e ) and W s (M. e ) 
intersect transversely at {y = 0}, with the transversality being 0(e^ n+2 ^ 2 ). Furthermore, 
the time the resultant trajectory spends in the region < x < n is given by 

T f = 0(2raln-). 

5.4 Completion of argument 

Using the notation of the previous section, near the invariant plane {x = 0} write 

W u (M t ) = {(x, y,z,e) : x = U x (y, z,e),y- + 77 < y < y+ - n} 
W s (M e ) = {(x, y,z,e) : x = S x (y, z, e), y_ + rj < y < y + - rf\. 
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By Theorem |5.19 it is known that if < p — //(e) = 0(e n ) for n > 3 and if the pair (di, cfo) 
satisfy the hypotheses of Lemma 5.17| , then there exists a (5(e), 0,0) G Co, x(e) = 0(e ra / 2 ), 
such that 

W x (0,0, e) =5 x (0,0,e) 

sgn((W* - <S£)(0,0,e)) = sgn(aidi + Q 2 d 2 ) 

(Wf - S x z )(0, 0, e) = 0(e(™+ 2 )/ 2 ). 

Since the manifolds intersect transversely at (x(e),0,0), there is an orbit, (p(t) , u(t) , <fi(t)) , 
which connects the two critical points (x*(e), 0, ±z*(e)) G M. e . Let this trajectory by 
translated so that 0(0), u(0), 0(0)) = (x(e),0,0) G C . 

Since the manifolds W"(A^ e ) and W s (Ai £ ) intersect transversely at {y = 0}, they 
intersect transversely everywhere along the orbit. The next goal is to understand the 
transversality as W u (M e ) intersects . Define 

6 = ( P '(o),u'(o),<p'(o)) T 

6 = (U*(x(e),0,l) T (5.37) 
£3 = (Sf(x(e),0,l) T 

It is clear that £ 2 • * G TVF M (7W e ) and £3 • i G TW 75 ^,) for all i, and that & • t G 
7W n (.M e )n7W s (.A/f e ) for all t. The three-form 

-fzyz will be used to gain an understanding 
as to how the flow carries these vectors. Using ( |5.29| ), the three-form P xyz satisfies Abel's 
formula 

P' — —37/ P 

1 xyz ~ OUJr xyzi 

which, using the fact that u = p'/p, has the solution 

P xyz (t) = ?^Lp xyz (0). (5.38) 

When applied to the tangent vectors £j defined in (]5 .37j ) , 

P xyz (0) = (A(0) + 0(e 2 ))(^ -5f)(0,0,e); (5.39) 
thus, after substituting into (|5.38|) one gets 

P xyz (t) = ^|y(A(0) + 0(e 2 )){W z -«S*)(0,0,e). (5.40) 

Now define T u such that p(T u ) = xq — v. As they intersect the section {x = xo — v} the 
manifolds W u {M e ) and W s {M t ) are given by the curves 

W u (M e )nB+ = {(y,z,e) : y=W(z,e)} 
W s (M e )DB+ = {(y,z,e) : y = Sy(z,e)}. 

In addition, there exists a z(e) such that 

In order to understand the transversality of the manifolds at {x = xq — v} it is necessary 
to compute (Wf — Sf )(z(e), e). Towards this end, set 

£2 = (0,^(z(e),e),l) T 
f 3 = (0„Sf(*( e ),e),lf. 
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It is clear that 



= p'(T v )(Uy-Sy)(z(e),e). 
Since p'(T u ) > is 0(u), upon substituting ( |5.40| ) into fl5.41| ) one gets 
(Uy-Sy)(~z(e),e) oc p 3 (0)A(0)(^' - Sf)(0, 0, e) 

= A: 2n+1 +0(e 2n+2 ), l^) 

where sgn (A) = — sgn (aidi + 02^2) (recall that A(0) < 0). The argument for the following 
lemma is now complete. 

Lemma 5.20 Let < p — p(e) = 0(e n ) for some n > 3. There exists a z(e), with 
\z(e) - <f>(T v )\ = 0(e n ~ 1 ), such that 

1. l(v(z(e),e)=SV(z{e),e) 

2. sgn — <SJf)(5:(e), e)) = —sgn (a%d\ + 02^2) 

3. (Uy-Sy)(~z(e),e) = 0(e^). 

Remark 5.21 To paraphrase, the manifold W u (Ai e ) intersects W s (M e ) transversely at 
the section {x = xq — u}, with the transversality being 0(e 2n+1 ) for < p, — //(e) = 0(e n ). 

6 Existence of solitary waves 

6.1 Existence of bright solitary waves 

A bright solitary wave is characterized by 

lim x(t) = 0; 

t— >±oo 

hence, in order to have such a wave it is necessary that W u (0) n P^ s (0) be nontrivial. Due 
to the symmetry described in Proposition 2.1, it is enough to show that W u (0) n C xo ■ t xo , 



and hence ^"(0)11^. 

Recall Lemma 13.11 , which states that C Xo ■ t xo is C 1 — 0(e~ C//<E ) close to W s (Ai e ) for 



points p £ C xo which take 0(1/ e) time to exit B xo . Under a time reversal, the flow on A4 e 
for z = O(e) satisfies 

z' = e 2 a + 0(e 3 ). 

Therefore, C xo -t xo will be close to W s (M. e ) in a region where z = 0(e) with sgn (z) = sgn (a). 

Now recall Corollary |5,4| , which states that when p = p(e) that <j>(T v ), the z-coordinate 
of W u (0)nB+ , is given by 

<f>(T v ) = -x (dt + ^x 2 d 2 + 0(v))e + 0(e 2 ). 

Thus, by the comments of the previous paragraph, if 

3 

(di + -xld 2 )(J < 0, 

then W u (0) n W S (M € ) is C 1 -0(e~ c / e ) close to the curve C Xo ■ t Xo (see Figures | and |). 
Using the facts that C xo ■ t XQ is C 1 -0(e~ c / e ) close to W s (7W e ) and that 1^(0) intersects 
W s (M t ) transversely yield the following lemma. 
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Lemma 6.1 Suppose that 

3 

(di + -x\d 2 )a < 

(see Figure^). There exists a /^(e) < n(e), with /x(e) — /U/j(e) = 0(e _C//<E ), smc/i that when 
[i = /ift(e), then W u (0) D W s (0) is nontrivial. 



Remark 6.2 The fact that //(e) — //^(e) > is a consequence of Lemmas 4J_ an d 4j_£- 

The solution described in the above lemma can be thought of as a 1-pulse solution. 
It is characterized by the fact W u (0) intersects the plane {y = 0} at exactly one point. 
An iV-pulse solution is characterized by both W u (0) n VF s (0) being nontrivial and W u (0) 
intersecting {y = 0} at N distinct points. Given the 1-pulse solution, it is natural to inquire 
as to the existence of iV-pulses. Fortunately, this question has been studied in the work of 
Kapitula and Maier-Paape [?]. In order to quote the results stated in that paper, the next 
lemma must first be proved. 

Lemma 6.3 Let < e <C 1 be fixed. Let Zh(n) represent the z-coordinate of W u (0) n {y = 
0}. Then 

1. z h (n h (e)) = 

2. J^ z Me)) + 0. 

Proof: The first part of the conclusion follows immediately from the fact that when /i = 

ft(e)T(0)n^^. 

For fixed e the set W u (0) n B+ Q yields a curve parameterized by fi, i.e., 

W u (0)nB+={(y,z,fx) :y = y u ^),z = z u ^)}. 



Due to Lemma 4.1 it is known that 

t/w > o, 

so that 

^-y u Me)) > 0, (6.1) 
dfi 

as Mh (e) =0( M (e)) = 0(e 2 ). 

Let p = (x - u,y u (fj, h (e)),z u (fi h (e))) G 1F U (0) n B+, and let T h > be such that 
p • Th £ C xo . Since = — e 2 c + 0(e 3 ) and z u {^h{e)) = O(e), it is necessarily true that 
Th = 0(1/ e). Therefore, as a consequence of the Exchange Lemma with Exponentially 
Small Error [?], the transversality described by ( |6.1D gets transferred into a transversality 
condition in the slow direction. Since the slow direction near M e is described by z, this 
means that 

—z h (n h (e)) / 0. ■ 

With the above lemma in hand, it is now possible to state Theorem 1.7 in Kapitula and 
Maier-Paape [?]. 
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Theorem 6.4 Let < e < 1, and suppose that 

3 

(di + -xld 2 )a < 0. 

For each N > 2 there exists a bi-infinite sequence {^} such that when fi = fj,^ there is an 



N -pulse solution to (2.7). The N -pulse is such that x(t) is even in t. Furthermore, 

\^-Me)\ = 0(e- c ^) 

as \k\ — > oo. 

Remark 6.5 The estimate on — ///i(e)| is not explicitly provided in Theorem 1.7 of 
Kapitula and Maier-Paape; however, it is implicit in the proof of that theorem. 

Remark 6.6 Although one can discuss the existence of N -pulses which are odd in t, this 
will not be done here. For a more complete description of the dynamical behavior for [i near 
/i/i(e), the interested reader should consult [?]. 

6.2 Existence of dark solitary waves 

A dark solitary wave is characterized by 

lim (x,y,z) • t G M € ; 

hence, in order to have such a wave it is necessary that both W u (M e ) n W s (M e ) ^ and 



there exist critical points in Ai £ - By Proposition ^4, the existence of the critical points is 
guaranteed if 

(di + 2xld 2 )a < 0. (6.2) 
By Lemma W u (M e ) n W s {M e ) ^ if 

0< M - M (e) = O(e"), n>3 
(di + xld 2 ){aid 1 + a 2 d 2 ) > 0, 



where, by Lemma 5.5, the coefficients oci satisfy 



Furthermore, the manifolds intersect transversely, with the transversality being 0(e 2n+l ) 
(Lemma 5.2C| ). 



An A^-circuit solution is a dark solitary wave whose trajectory passes near {x = 0} N 



times. If (6.2) and ( |6.3j ) are satisfied, then for < e <C 1 there is a 1-circuit solution. The 
goal is to show that there exist iV-circuits for 2 < iV < N(e), where N(e) — > oo as e — > 0. 
This will be accomplished through the use of the Exchange Lemma with Exponentially 
Small Error (ELESE) (Jones et al [?]). 

In order for the ELESE to apply, it must first be shown that W u (Ai € ) D W s (M e ) 
transversely. This task was accomplished in the series of lemmas leading to Lemma |5.20| . 
The ELESE then applies to all those points in 

C+ = W u (M e )nB+ 
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which take 0(l/e) time to exit B xo . 
Set 

C+ = W s (M e )B+. 
The z-coordinate of n C+, say z, is given by 

3 

z = -x {di + -x\d 2 + 0(v))e + 0(e 2 ), 

i.e., \z — <j>(T v )\ = 0{e 2 ). This is due to the assumption that \i — //(e) = 0(e n ) for some 
n > 3. Set 

c u - = r(M £ )ns; , c.- = W(A< e )niJ-. 

Due to the symmetry inherent in the ODE, not only is C~ f~)C~ 7^ 0, the z-coordinate of 
the intersection is given by —z. Now, recall that for \z\ = O(e) the flow on A4 e is given by 

z' = -e 2 a + 0(e 3 ) 

(equation 3.39|) . Thus, if 



(di + \x 2 Q d 2 )a < (6.4) 

the flow on M e is from z to —z; otherwise, it is in the opposite direction. 

Assume that ( |6.4D holds. Since z' = 0(e 2 ), it takes a time of 0(l/e) for trajectories 
to traverse from z to —z. Let ■ represent the intersection of the curve with the 
section B~ Q as it is carried by the flow generated by the ODE (^). By the ELESE the 
curve C+ • T+ will be C l -0{e~ c / e ) close to C~ at —z. Since C u D C s transversely at —z, 
with the order of transversality being 0(e 2n+1 ), this implies that • T+ n C~ transversely 
exponentially close to —z. In other words, upon passing through B XQ the manifold W u (A4 e ) 
becomes oriented so that it once again intersects W s (A4 e ) transversely. Given the existence 
of a 1-circuit solution, the existence of a 2-circuit has now been proven. 

The above argument can be repeated to show the existence of iV-circuit solutions for 
each N > 2. Given an e > 0, the maximal N, say N(e), will be such that N(e) — > 00 
as e — > 0. This is due to the fact that upon passage through B xo , the curve • T+ is 
C 1 -0(e~ c / e ) close to C~ at —z. The following theorem has now been proven. 



Theorem 6.7 Suppose that < f « 1, and suppose that (6J.), (6.6 ), and l \6.4) hold (see 
Figure Then there exists an N(e) > 1, with N(e) — > 00 as e — ► 0, such that N -circuit 
solutions exist for 1 < N < N(e). 
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Figure 1: Flow on {z = 0} when e = 
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Figure 2: The phase space geometry 
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Figure 3: Existence of dark solitary waves 
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Figure 4: Existence of bright solitary waves 
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Figure 5: Projection of flow on {y = 0} 
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Figure 6: Parameter regime for bright solitary waves 
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Figure 7: Parameter regime for dark solitary waves 
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